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SUMMARY 

A method 1b presented for determining the load distribution due 
to steady roll and pitch for thin flat-plate wings whose plan form 
is arbitrary except that a part of the leading edge must be supersonic 

For wings with straight supersonic leading edges, the load 
distributions due to angle of attack, steady roll, and steady pitch 
are explicitly evaluated and are computed for a family of wings whose 
plan form includes most types of flow region commonly encountered. 
These computations showed that negative lift existed toward the rear 
of pointed wings whose aspect ratio was small. In steady roll, 
negative loading occurred In regions influenced by the edge of the 
plan form at the opposite side of the roll axis. When the pitch 
axis was located near the semichord position, the load gradient for 
steady pitch was approximately in the chordwise direction, except 
in regions Influenced by subsonic trailing edges. High positive 
loading occurred toward the front of the wing and high negative 
loading toward the rear. 


INTRODUCTION 

A method is presented in reference 1 for determining the pressure 
distribution over thin wings at supersonic speeds. The method is 
based an an integration of the local source strength (which is 
proportional to the local slope of the wing surface) over the regions 
of the disturbed flow field that lie within the forward Mach cone 
from a point on the wing surface. Reference 1 shows that the 
contributions to the pressure coefficient of the disturbed fields 
off the surface of the wing may be replaced by equivalent contribu- 
tions of parts of the wing surface, and that some of the surface 
integrals that are involved in the determination of pressure coeffi- 
cient may be reduced to line Integrals. 

The reduction of surface integrals to line integrals is feasible 
for all regions of a flat-plate wing except those influenced by 
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interacting disturbed flow fields off the wing plan boundary. For 
the wing shown in figure 1, for example, all regions may be treated 
by the methods of reference 1 except the small shaded regions at the 
rear. The restriction that a portion of the leading, edge must be 
swept ahead of the Mach lines from the foremost point of the wing 
thus guarantees that some portions of the wing will be subject to 
the methods of reference 1. For regions influenced by interacting 
disturbed flow fields, more elaborate methods are required, such as 
those used for delta wings in reference 2. 

The load distributions due to roll and pitch have been determined 
for seme plan forms with straight edges in references 3 and 4. The 
methods of reference 1 are applied herein to the determination of 
these load distributions for more general classes of plan form, whose 
edges may be curved . For a family of wings of the type shown in 
figure 1, but having, for convenience, straight supersonic leading 
edges, the load distributions due to angle of attach, steady roll, 
and steady pitch were computed. This type of wing was chosen 
because it contains most types of flow field commonly encountered. 

This analysis was completed at the MCA Cleveland laboratory 
during January 1948. 


SYMBOLS 

The following symbols, some of which are illustrated in fig- 
ure 2 to 6, are used throughout this report: 

A, B, Ci, 

_ „ substitution terms 

c 2 > • • •> U 10 



integration limits 
pressure coefficient 

slope of straight leading edge in (u,v) coordinate 
system 

Mach number 

steady rate of roll, radians per second 
steady rate of pitch, radians per second 



C/fi 
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area or area integration. 

area on wing surface 

free -stream velocity (in x-directlon) 

variables of integration in oblique coordinates 

oblique coordinates of point on wing surface 

oblique coordinates of wing vertex 


functions defining form of right supersonic leading 
edge 

functions defining form of left supersonic leading 
edge 

functions defining form of right subsonic leading 
and trailing edges 


functions defining form of left subsonic leading 
and trailing edges 

component of perturbation velocity in z-dlrectlon 
(positive outward from z ■ 0 plane) 


y 

n 

a 

3 

a T> 


Cartesian coordinates of point on wing surface 
variables of integration in Cartesian coordinates 
coordinate of roll axis 
coordinate of pitch axis 
angle of attack 



w/u for top and bottom wing surface, respectively, 
in y a constant plane 
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ANALYSIS 

A plan form that contains most types of wing region commonly 
encountered is shown in figure 1. The downstream Mach lines from 
the vertex and from the junctures of the subsonic and supersonic 
sections of the leading edges divide the plan form into nine types of 
region that differ in the type and number of wing edges that affect 
the pressure distribution. Regions I and II are influenced only by 
supersonic leading edges, whereas regions III, IV, and V are affected 
also by one subsonic leading or trailing edge. In regions VI, VII, 
and VIII, the subsonic edges of both sides of the wing affect the 
flow. The shaded areas represent regions that are affected by 
Interacting perturbed fields off the wing plan form and are not 
easily treated by the methods of reference 1. 

The essential equations required to determine the load distribu- 
tion for the type of wing shown in figure 1 (without the shaded 
regions) may be obtained from a consideration of regions of types III, 
IV, and V. For these types of region (influenced by supersonic 
leading edges and only one subsonic leading or trailing edge), an 
expression was derived in reference 1 for the pressure coefficient 
at a point (x,y) when that point is not influenced by vorticity off 
the wing plan-form boundary. This expression is (see fig. 2): 
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where ia the derivative of the equation defining the plea 

form, of the subsoalc leading and trailing edges evaluated at the 
intersection of the right forward Mach line from (x,y) and the iring 
boundary (point d). Equation (l) was shown in reference 1 to contain 
an additional tern of a form similar to the line integral along the 
Mach line bd. This additional line integral is related to the 
vorticity behind the trailing edge and is therefore zero for leading 
edges. For trailing edges, the nature of the line integral depends 
on the conditions Imposed. In particular, if the Eutta-Joukowski 
condition is imposed, the additional line integral must be such that 
it exactly cancels the integral along the Mach line bd. Henoe, 
equation (1), -without the last integral, can be used to determine the 
pressure coefficient in regions influenced by vorticity off the plan 
form provided that the Eutta-Joukowski condition is imposed at 
subsonic trailing edges. The imposition of this condition, although 
arbitrary, is conventional and will be assumed in the rest of the 
analysis. 


For a flat-plate wing, the effective local slopes O b and Otj> 
tor determining the loading due to angle of attack, steady roll, and 
steady pitch, are 


a 


T 


-a 


a B » a 


^ § (tj - n 0 ) 

a B - § <n - no) 




T 


a B 


1 


u 


ce- 1 0 ) 


*q) 


( 2 ) 


(2a) 


(2b) 


where p and q are the rates of roll and pitch, respectively , in 
radians per second, and i)q and £q are the distances from the 
origin of coordinates to the roll axis and the pitch axis, 
respectively. From equations (l) to (2b) the first, second, and 
fourth integrals of equation (1) are seen to vanish for the pressure 
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coefficient due to angle of attack and roll, whereas the second and 
fourth integrals vanish for the pressure coefficient due to pitch. 

It is convenient to convert the required integrals to an oblique 
coordinate system whose axes are the rearward Mach lines from the 
origin. The origin of the coordinates is taken at the juncture of 
the supersonic and subsonic parts of the leading edge. The conversion 
equations are (fig. 3) 


u 88 ~ Pn) 

i m |j(v + U) 

Uw “ 2p( x 

x = R< v w + «w) 


v - + Pq) 

Tl - |(T - U) 

v w = ||< x + py> 

y * r^w - 


(3) 


The elementary area in this coordinate system is ^ dudv. The 

coordinates (v^, v w ) or (x,y) are used to represent the point 
on the wing for which the pressure coefficient is desired, whereas 
(u,v) or (|,tj) represent variables of integration. In the oblique 
coordinate system, the equations for the load distributions due to 
angle of attack, roll, and pitch become, respectively, (fig. 4) 





b 

2 dv-du 

* ^(i^-u) (v w -v) 

-'a 


(4) 


pn 0 i 



(v-U3-Mno)dv 

/\Jv^ 


2 

jtMq 0 


rb 


(v-u-Mtjq) (dv-du) 

l||(uw- u )(v w -v) 


(4a) 


a 
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W (y+u 2 

rfeo)»T 

/f«w- u 3 


v v 


gp (™-fQ(dy-du) 



(4b) 


where u* and Its derlyatlyee are eyaluated far y * v w . Equa- 
tions (4). (4a), and (4b) can be integrated for arbitrary forms 

of u 5 (y v ), If u and v are linearly related along the line ab. 

In other words. If the supersonic part of the leading edge Is a 
straight line, explicit expressions for the load distribution are 
readily obtainable for arbitrary forms of the subsonic leading and 
trailing edges. If the supersonlo leading edge Is not a straight 
line, the load distributions due to angle of attack and roll may 
still be obtained by means of a simple graphical Integration along 
the required part ab of the leading edge. The treatment of such 
wings to determine the lift distribution is described In detail In 
reference 5. This method may be readily extended to determine the 
line Integrals along ab for load distribution due to roll and pitch. 
The area Integral required for the pitch loading (equation (4b)), 
although somewhat more difficult to evaluate than the line Integrals, 
is also subject to stripwlse , graphical Integration methods. 

The Integrations from b to v w in equations (4) to (4b) 

(along the Mach line u = u 3 (v w )) are Independent of the form of 
the wing boundary and hence may always be explicitly Integrated. 

If the Kutta-Joukowskl condition is Imposed at the trailing edge, 
these integrals need be evaluated only when the right forward Mach 
line from (u w ,v w ) Intersects the plan-foam boundary at a subsonic 
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leading edge ( 0 < < 1 1. If the right forward Mach line Inter- 

' / / dug \ 

sects the plan-fora houndary at a subsonic trailing edge lJ j 

the integrals along u ■ UjCv^) vanish. 


The analysis may now be extended to regions of the wing 
Influenced by all its subsonic leading and trailing edges. In general, 
suoh regions may be of two types (figs. 5(a) and 5(b)). In both types 
of region, the left forward Mach line from (u w ,v ir ) intersects a 
subsonic leading or trailing edge at the opposite side of the wing. 

In figure 5 (a), (u^v^) lies In a region of type 71 (fig. 1 ), tar 

which the reflections of the forward Mach lines at the wing boundaries 
do not intersect on the wing surface. In figure 5(b), («*,▼*) lies 

in a region of type VIII, for which the reflections of the forward 
Mach lines cross on the wing surface. For both types of region, the 
methods of reference 1 indicate that an additional line Integral 
along the Mach line v « T&tu,,) must be added to equation (4). These 
integrals are the same as those along u * U 3 (▼,,), except that y 
replaces u, and V 4 replaces u 3 . The sense of the integration Is 
again from the supersonic leading edge V]_ or Yg to the subsonic 
edge Y4. The integration along y * V4, like that along u « u 3 , 
vanishes (for solutions that satisfy the Khtta- J oukowski condition) 

If the forward Mach line from (u w ,y w ) Intersects a subsonic trailing 
edge (as it does, for example, in fig. 5(b)). 


Along the supersonic leading edge, the sense of the Integration 

is from a to b for both types of region. The values of the line 

Integrals along ab are thus of opposite sign for figures 5(a) and 5(b). 

The area integrations for the pitch loading (equation (4c)) extend 
over the shaded areas of figures 5(a) and 5(b). For the type of 
region shown in figure 5(b), the area integration consists of two 
parts. The integration over the downstream area is Independent of 
the fora of the supersonic leading edge. The integration over the 
upstream area depends on the contour of the supersonic leading edge. 

This integration is subtracted from the integration for the lower area, 
because two area cancelations are involved, one for each of the regions 
off the two subsonic leading edges (reference l). 

For regions of type VII (fig. l), both right and left reflected 
Mach lines Intersect the plan-form boundary along the same supersonic 
leading edge. A point in this region is like a point in region VIII 
if the Mach lines cross on the surface and is like a point in 
region 71 if the Mach lines do not cross. 



MCA TN Ho. 1689 


9 


For regions of types I and II, both forward Mach lines intersect 
the plan-fora boundary at supersonic leading edges. Hence, all of 
the line integrals except that along the supersonic leading edge 
vanish, and the limits for this integral are the points of intersec- 
tion of the forward Mach lines with the leading edge. The area 
integration for the pitch loading extends over the region bounded 
by the forward Mach lines and the leading edge. 


APPLICATION TO WINGS WITH STRAIGHT 
SUPERSONIC LEADING EDGES 


When the wing is symmetrical and a section of the leading edge 
on both sides of the line of symmetry is straight and swept ahead 
of the Mach lines from, the vertex, the equations for the supersonic 


leading edges are U]_ 


kv and U2 



The origin of 


coordinates is taken at the juncture of the supersonic and subsonic 
leading edges (fig. 6). The axis of symmetry is identified with the 
roll axis. For a general region on such a wing, line Integrals are 
required along the Mach lines reflected from the subsonic leading 
edges and along the two sections of the supersonic leading edge. 

The limits for these line integrals , as well as those for the area 
integration required for pitch loading, vary with the type of region 
considered. The presentation is simplified if the expressions for 
a general region are first evaluated and the appropriate limits for 
each region are then indicated. 


In the following expressions for the load distributions due to 
angle of attack, steady roll, and steady pitch, respectively, the 
first integral is the line integral along v * v 4 (u) , the second is 
the integral along u = u 3 (v), the third is along u = u^(v) , and 
the fourth is along u * ugtv) : 


a 


2 




1 ^ 


a 




dy 



1+k 

rc , 

— (l+k) 

f b du 

k 

b» j 

^(u^-u) [v w +ku-M»iQ ( l~fc)J 


( 5 ) 
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j * where 



is the integration over the required areas for each region. 

From these integrals, the following explicit expressions are 
obtained for the load distributions due to angle of attach, steady 
roll, and steady pitch, respectively: 



( 6 ) 
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973 



973 



(1-k 2 ) m 0 - 8 *££ 0 ] - 



14 


MCA Ttl Ho. 1689 


The appropriate limits for each of the eight regions shown in 
figure 6 are given in the following table: 


Limit 

Wing regions (fig. 6) 

I 

II 

hi 

IV 

V 

r — — — — - — • — 

VI 

VII 

VIII 

a 

Sf 

Sr 

Sr 

Sr 

Sr 

v 4 
" k 

v 4 
" k 

-kv 4 

a 1 

Sr 

Sr 

Sr 

Sr 


Tr (1 - k) + 

v 4 

k 

>»* 

v 4 

k 

u d 

b 

Sr 

Si 

Sr 

u d 

u 3 

u d 

u 3 

u 3 

b' 

Sr 

u d 

Sr 

u d 

u 3 

u d 

u 3 

-kv 4 

c 

-kv w 

-kSr. 

u 3 

u 3 

u 3 

u 3 

u 3 

u d 

c 1 

Sr 

"Sr 

-u 3 /k 

< 

to 

a 

I 

-ku 3 + 
Mn 0 ( i -k) 

-u 3 /k 

-ku 3 + 
Mt| 0 (l-k) 

-ku 3 + 
Mrj 0 (l-k) 


An examination of equations (6) to (6b) shows that the arc-tangent 
factors are the same for each type of load distribution. The square- 
root terms are also the same for the roll and pitch loadings. 
Application of the foregoing table to equations (6) to (6b) shows 
that many terns vanish for some of the wing regions because the upper 
£hd lower limits are identical. The location of these limits for 
regions of types VI and VIII are shown in figure 5. 

The expressions for S in each region are: 


07 ^ 
















Area integration, S 


11- 

Cl J A 

+ bc 2 



2C 3 

+ ac 4 




2C 3 

+ a(c 4 

- Ci) + BCg 


2C 5 

+ bc 6 




2(C 7 

+ C 3 ) 

+ A(c 4 - 

c l) 

+ B(C, 

2(C ? 

+ c 5 ) 

+ b(c 6 - 

c 8> 


2(C 5 

+ c io' 

) + A(c ± 

- Cg) 

+ B(i 


(VU-VI3) [v w +ku 3 -Mn 0 (!-*)] 


Y w 4ku 3 -Mn 0 ( 1-k ) 


C£ ■ tan 

C 3 = \f( u w -u 3 )'(t v +^) 


ku v +V4»Mrio (l-k) 


t w“ t 4 
1 u w .+krv , 4 


C 4 - — ' \| kv..+u 


u^+k^) (v v -v 4 













16 


MCA TIT No. 1689 


Equations (6) to (6b) have been used to calculate the load 
distributions due to lift, roll, and pitch for the wing shown in 
figure 7, For computation, use of the coordinates 

and v w ^iJtiq| rather than was convenient . .These coordinates 

make equations 6 to 6(b) nondimensional . 

Ph 0 | 

For the wing shown in figure 7, the ratio ~j was taken 


equal to 1.0. The value of k is l/3. 
subsonic edges were assumed to be 


The equations for the 


u 3 


M|-l 0 | 


l Chol) 


T 4 


«kl 


1 + 4 


ter ■ ‘)‘ 


(7) 


In order to satisfy the Kutta-Joukowski condition, the integral 
along u * uj is zero for v w /(m|t[q|) > 12 “^ 2 ; and the integral 

along v « v 4 is zero for u w / (M |tj 0 | ) > 1 + 12" 1 / 2 , because 
du 3 dv. 

and are then greater than unity. 

The contour of the wing is represented in figure 7 for a Mach 
number of *f2, although in the (u,v) cool'd inate system the plot 
represents a series of wings whose spatial contours vary with Mach 
number according to equations (7) and the value k = l/3. Hence 
the load distributions calculated for this wing apply directly to, 
all wings of the series defined by k = l/3 and equations (7). The 
load distributions for a considerable variety of plan forms, at a 
given Mach number, can be obtained from the load distributions 
calculated for the wing of figure 7 by terminating the wing with any 
form of supersonic trailing edge. The load distributions for the 
remaining regions of the wing are unaffected by such changes. 

The effect of altering the location of the roll or pitoh axes 
can be determined with the aid of the superposition principle. If 
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the roll and pitch axes are shifted to and i lf respectively, 
then equations (2a) and (2h) "become, for roll, 

a T “ " “ £(t) 0 -i}i) 

a B = + u^O -T il) (8) 

and, for pitch, 

°T = ' 


a B '§U-?o) «*) 

Because ^‘( T lo“ T il) u^O"^l^ are constants, the contributions 

to the load distributions due to these terms are exactly equivalent 
to the load distributions due to the corresponding angles of attack. 
Thus, if the load distributions are computed for the axes t}Q 

and Iq (or for some relation such as Pj^o |/^0 “ as assumed 

for the wing of fig. 7), the load distributions for roll or pitch 
about any other axes are simply the load distributions for the 
axes t|q and i Q plus the lift distributions for the angles of 
attack: 

a ■ * uW 


and 


a = ± §( ! 0 “ 

For a family of wings whose contour is represented by the 
value k = l/3, and equation (7), the load distributions due to 
angle of attack, steady roll, and steady pitch are shown in fig- 
ures 8, 9, and 10, respectively. The dashed lines in these figures 
are Mach lines that separate the various regions indicated in 

figure 7. The additional Mach lines at Uy/(M|ilQ|) = 1.29 and 
at v w / ( m|tjo | ) * 0.29 separate the regions influenced by the subsonic 
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trailing edges from those influenced only by leading edges. The pres- 
sure coefficient is zero along the subsonic trailing edge for each 
type of loading, because the Khtta-Joukowski condition was imposed 
in each case. Along each subsonic leading edge the pressure coeffi- 
cient is infinite. 

The lift distribution shows that positive lift exists on all 
parts of the wing surface except the extreme rearward tip (region TUI 
and part of region VII). This negative region is a result of the 
upwash over the subsonic edges. The lift decreases rapidly toward 
the subsonic edges. The region of the wing having the highest lift 
is that bounded by the leading edges and by the Mach lines from the 
beginning of the trailing edges; 

For steady roll (fig. 9), the loading is negative in the outboard 
part of region IV and in nearly all of regions V, VI, and VII. In 
region VIII, the loading again becomes positive. The large negative 
region results because the greatest contribution to the loading 
proceeds from the leading edge on the opposite side of the wing, 
where the vertical component of the perturbation velocity is of 
opposite sign. 

The load gradient for steady pitch (fig. 10) is primarily in 
the chordwlse direction, except in the regions influenced largely by 
subsonic trailing edges. High positive loading oocurs toward the 
front of the wing and high negative loading toward the rear. The 
loading becomes negative ahead of the pitch axis because the contribu- 
tion due to the area integration and the contribution due to the line 
integrals are of opposite sign ahead of the pitch axis (|q > |). The 
loading therefore changes sign when the contribution of the Integrals 
over the area included in the forward Mach cone is sufficiently 
large to overbalance the contribution due to the line integrals. 


SUMMARY OF THEORY AND RESULTS 

A method has been presented for determining the load distribution 
due to steady roll and steady pitch on thin wings whose plan form is 
arbitrary except that a part of the leading edge must be supersonic. 
When the supersonic part of the leading edge is a straight line, 
these load distributions can be explicitly evaluated for all regions 
of the wing except those influenced by interacting flow fields off 
the plan form. 

For a particular family of wings having a plan form that 
includes most types of flow field commonly encountered, the load 
distributions due to angle of attack, steady roll, and steady pitch 
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wore computed. The lift distribution for this family of wings 
showed that negative lift may exist toward the rear of pointed 
wings if the aspect ratio is small. The highest lift occurred in 
regions affected only by leading edges. In steady roll, negative 
loading occurred in regions influenced by the edge of the plan form 
at the opposite side of the roll axis. At the extreme rear of a 
low-aspect-ratio wing, the loading again became positive. With the 
pitch axis located near the semichord position, the load gradient 
for steady pitch was primarily in the chordwise direction except 
in regions Influenced by subsonic trailing edges. High positive 
loading occurred toward the front of the wing and high negative 
loading toward the rear. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 15, 1948. 
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Figure 8. - Lift distribution for series of wings defined by k = l/3 and equations (7). 
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